PHOTOGRAPH  THIS  SHEET 


i**i 


& 


LEVEL 


m 

INVENTORY 


FISD-TK-77-S?  ' 

DOCUMENT  IDENTIFICATION  ^  7 7 


g  Ed'tnondSj  iJ  B*; 
o  I - 


f?  u*77 


DISTRIBUTION  STATEMENT  A  ! 

Approved  for  public  release;  ! 

Distribution  Unlimited  i 

I _ - : _ i 

DISTRIBUTION  STATEMENT 


ACCESSION  FOR 
NTIS  GRAAI 

DTIC  TAB 

UNANNOUNCED 
JUSTIFICATION 


BY _ 

DISTRIBUTION/ 

AVAILABILITY  CODES _ 

DIST  I  AVAIL  AND/OR  SPECIAL 


DTIC 

ELECTE 
APR  261963 


D 


DATE  ACCESSIONED 


copy 

(INSPECTED  , 

x.  2  y 


DISTRIBUTION  STAMP 


83  04  26  051 

DATE  RECEIVED  IN  DTIC 

PHOTOGRAPH  THIS  SHEET  AND  RETURN  TO  DTIC-DDA-2 


FORM  Jk 

DTIC  OCT  79  70A 


DOCUMENT  PROCESSING  SHEET 


8  I Z  l  Z  TV  Q\f 


ASD-TR-77-59 


MULTIVARIABLE  CONTROL  FOR 
A  VARIABLE  AREA  TURBINE  ENGINE 


PERFORMANCE  /  STABILITY  DIVISION 
DIRECTORATE  OF  PROPULSION 
PROPULSION  DEPOT  ATE 


AUGUST  1977 


TECHNICAL  REPORT  ASD-TR-77-59 

Final  Report  for  Period  January  1976  —  March  1977 


Approved  for  public  release;  distribution  unlimited. 


PROPULSION  DEPUTATE 
AERONAUTICAL  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  ay  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office  (10)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NTiS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


David  B.  Edmunds 
Aerospace  Engineer 
Performance/Stability  Division 
ASD/YZEA 


FOR  THE  COMMANDER 


AJames  L.  Petti  g'rew,  It  Col,  USAF 
Chier,  °erformance/Stabi 1 i ty  Division 
Directorate  of  Engineering  and  Test 
Propulsion  Deputate 


Copies  of  this  report  should  not  be  returned  unless  return  is  required 
by  security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document. 

AIR  FORCE/56780/10  February  1978—  150 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

ASD-TR-77-59 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (md  Subtitle) 

A  MULTIVARIABLE  CONTROL  FOR  A  VARIABLE  AREA 

TURBINE  ENGINE 

S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Technical  Report 

January  76  -  March  77 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf*.) 

David  B.  Edmunds 

8.  CONTRACT  OR  GRANT  NUMRF.Rf» 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Aeronautical  Systems  Division  (AFSC) 
Performance/Stability  Division  (YZEA) 
Wriqht-Patterson  AFB,  Ohio  45433 

to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 

-328A 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Directorate  of  Engineering 

Deputy  for  Propulsion 

Wriqht-Patterson  AFB,  Ohio  45433 

12.  REPORT  RATE 

August  1977 

13.  NUMBER  OF  PAGES 

90 

14.  MONITORING  AGENCY  NAME  A  AODRESSfJ/  different  from  Controlling  Ottlca) 

15.  SECURITY  CLASS,  (o*  thte  report) 

UNCLASSIFIED 

15*.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thte  Reoort) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (o i  the  eb, tract  entered  In  Block  30,  II  dl Iterant  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  ( Continue  on  rereree  aide  if  neeeeeary  and  identify  by  block  number) 

Controls,  Multivariable  control.  Variable  cycle  engine,  Turbine  engine 

20.  ABSTRACT  (Continue  on  raver aa  aide  it  neeeeeary  and  identify  by  block  number) 

A  multivariable  control  has  been  formulated  for  a  turbofan  nonaugmented  engine 
with  variable  area  turbines  operating  at  sea  level  static  conditions.  The  non¬ 
linear  dynamic  model  was  represented  by  two  fourth-order  linear  models.  An  in¬ 
teractive  matrix  Ricatti  optimization  program  was  developed  which  automatically 
calculated  a  feedback  gain  matrix  and  displays  on-line  the  transient  response  of 
selected  key  parameters.  The  multivariable  controller  implemented  on  the  non¬ 
linear  model  required  modifications  to  accommodate  both  large  and  small  transient; 
gains  transitioning  and  control  saturation.  Transient  results  from  the  controller 

DD  t  j2nM73  1473  EDITION  OF  1  NOV 65  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dmtm  Entered) 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PACEfWnn  JP«t«  Entottd) 

are  discussed  including  the  effect  of  component  deterioration,  sensor  and  control 
positioning  errors,  and  installation  factors. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEflWi»n  Dttt  Enfnd) 


TABLE  OF  CONTENTS 


SECTION 

I.  INTRODUCTION  AND  BACKGROUND 

II.  VARIABLE  CYCLE  ENGINE  NONLINEAR  MODEL 

1.  Variable  Cycle  Engine  Concept 

2.  Nonlinear  Dynamic  Model 

3.  Variable  Cycle  Engine  Model 

III.  CONTROL  DESIGN  METHODOLOGY 

1.  Linearization  Technique 

2.  Matrix  Reduction 

3.  Control  Optimization  Technique  for  Linearized  Models 

4.  Interactive  Rlcatti  Equation  Solver 

5.  Optimal  Control  Law  Implementation  on  Nonlinear  Dynamic 
Simulation 

IV.  NONLINEAR  MODEL  APPLICATION  RESULTS 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

REFERENCES 

APPENDIX  A 


LIS  OF  ILLUSTRATIONS 


FIGURE  TITLE  PAGE 

1  Variable  Cycle  Engine  Schematic  .  5 

2  Original  FlOO  Engine  Schematic  and  Component  Performance 

Characteristics  .  7 

3  Nonl inear  Model  Calculation  Flowpath  .  9 

4  Variable  Turbine  Efficiency  and  Flow  Area  Factors . 14 

5  Steady  State  Control  Schedules . ‘.16 

6  Steady  State  Output  and  State  Values.  ...  .  18 

7  Comparison  of  Linear  Model  with  Nonlinear  Model  for  Small 

Step  Change  in  Fuel  Flow  .  . . 34 

8  Transient  Response  for  Various  Performance  Indexes . 46 

9  Schematic  of  Regulator  Control . 53 

10  Feedback  Gain  Transitioning  Logic  .  55 

11  Transient  Response  for  VCE  Nonlinear  Model 

(PLA  =  20  to  PLA  »  83) . 62 

12  Transient  Response  for  VCE  Nonlinear  Model 

(PLA  »  83  to  PLA  »  20) . 65 

13  Transient  Response  for  VCE  Nonlinear  Model 

(PLA  *  73  to  PLA  »  83) . 68 

14  RCVV  Actuator  &*ror  Effect  on  VCE  Nonlinear  Model  .....  76 

15  Turbine  Deterioration  Effect  on  VCE  Nonlinear  Model  ....  77 

16  RCVV  Actuator  Error  Effect  on  Original  FlOO  Nonlinear  Model  78 

17  Turbine  Deterioration  Effect  on  the  Original  FlOO 

Nonlinear  Model . 79 

18  Reduced  Feedback  Gain  Matrix  Affect  on  VCE  Nonlinear  Model.  80 


LIST  OF  TABLES 


TABLE  TITLE  PAGE 

1.  State,  Input,  and  Output  Vectors  for  Original  F100  Model  .  .  10 

2.  State,  Input,  and  Output  Vectors  for  VCE  Model  .  21 

3«  1 8th  Order  Linearized  Matrices  for  PLA  =  60  . 

4»  Eigenvalues  of  A  Matrix  for  PLA  =  60 . 30 

5.  State  Vectors  for  6th  and  4th  Order  Models . 32 

6.  6th  and  4th  Order  Linearized  Matrices  for  PLA  =  60  .  33 

7.  Redefined  State,  Input, and  Output  Vectors . 39 

8.  Performance  Index  (PI)  Weightings . 44 

9.  Feedback  Gain  Matrices . 58 

10.  Output  Parameter  Changes  for  Oifferent  Implanted  Faults  - 

VCE  and  Original  F100  Nonlinear  Mode'1 . 72 


v 


NOMENCLATURE. 


Gas  Turbine  Parameters 

TT25H  ~  gas  total  temperature  downstream  of  inside  diameter  fan  exit 
guide  vanes. 

TT25c  ••  total  temperature  of  gas  in  duct. 

TT3  -  total  temperature  of  gas  at  compressor  discharge. 

TT4  -  total  temperature  of  gas  at  burner  exit. 

TT^M  -  gas  temperature  at  inlet  to  high  pressure  turbine, 
rr^hi  -  high  turbine  gas  flow  total  temperature  after  energy  loss 
to  vanes  and  blades. 

HVlo  -  high  turbine  gas  flow  total  temperature  after  energy  loss  to 
seals  and  discs. 

TT45hi  -  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
vanes  and  blades. 

TT^51o  -  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
seals  and  discs. 

TT5  "  low  turbine  gas  flow  exit  total  temperature. 

TT6c  -  duct  exit  gas  flow  total  temperature. 

TT7  “  augmentor  exit  gas  flow  total  temperature. 

PT25c  -  duct  entrance  total  pressure. 

PT3  -  compressor  discharge  total  pressure. 

PT*»5  ~  inter-turbine  volume  total  pressure. 

PT7  -  afterburner  volume  total  pressure. 

PT6c  -  fan  duct  volume  total  pressure. 

N1  -  fan  physical  speed. 

N2  -  compressor  physical  speed. 
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NOMENCLATURE  (cont.) 


AJC  -  core  exhaust  nozzle  area. 

CIVV  -  fan  inlet  guide  vanes. 

RCVV  -  compressor  stator  vanes. 

WF  -  main  burner  fuel  flow. 

AJD  -  fan  duct  exhaust  nozzle  area. 

HTVPOS  -  high  turbine  nozzle  position. 
FTVPOS  -  fan  turbine  nozzle  position. 

FN  -  net  thrust.. 

V/FAN  -  airflow  through  fan. 

SMAF  -  fan  surge  margin. 

SMHC  -  compressor  surge  margin. 

CMVT  -  Constant  Match  Varying  Temperature. 
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SECTION  I 


INTRODUCTION  AND  BACKGROUND 

The  majority  of  present  day  control  systems  for  gas  turbine 
engines  have  been  derived  using  classical  control  techniques.  These 
techniques  provided  acceptable  control  systems  but  required  increas¬ 
ingly  greater  designer  skill  and  time  to  exploit  the  full  performance 
capabilities  of  the  more  recent  variable  geometry  afterburning  tur¬ 
bofan  engines  (F100-PW-100,  GE-F101).  Proposed  advanced  gas  turbine 
engine  cycles  will  have  many  variable  geometry  components  including 

fan  and  compressor  vane  angles,  high  and  low  turbine  flow  areas,  and 
fan  and  core  exhaust  nozzle  areas,  in  addition  to  fuel  flow.  This 
presents  the  control  designer  with  a  complex  control  system  involv¬ 
ing  multiple  inputs  and  multiple  outputs. 

Although  control  modes  for  such  cycles  have  been  derived  using 
classical  1  inear  analysis  techniques,  Beattie  (Reference  1),  the  appli¬ 
cation  of  modern  control  techniques  and  optimization  procedures  offers 
a  systematic  approach  to  handle  these  more  complex  multivariable  problems. 
For  this  reason  interest  in  the  application  of  modern  control  theory  to 
the  design  and  analysis  of  jet  engine  controls  has  increased  in  recent 
years.  Bowles  (Reference  2)  and  Merrill  (Reference  3)  studied  the 
application  of  modern  control  theory  to  a  single  spool  turbojet  engine. 

Its  application  to  a  two  spool  turbofan  engine  has  been  studied  by 
Michael  and  Farrar  (References  k  and  5)  and  Weinberg  (Reference  6). 

Finally  Beattie  (Reference  7)  has  applied  multivariable  optimal  control 

techniques  to  a  variable  cycle  engine  (VCE) . 

These  studies  have  demonstrated  that  the  control  laws  derived 

from  the  optimal  control  theory  can  be  adapted  as  a  control  system 

on  dynamic  nonlinear  turbine  engine  models.  In  all  cases  the 


adaptive  control  system  was  suboptlmai  and  equal  or  better  than  the 
conventionally  designed  control  in  meeting  the  basic  requirements 
of  steady  state  performance  and  transient  response.  The  significant 
conclusion  from  these  studies  was  that  control  systems  derived  using 
modern  control  techniques  can  reduce  the  design  time,  trial  and  error 
and  intuition  required  relative  to  classical  control  techniques.  In 
addition  all  significant  interactions  between  inputs  were  accounted 
for  in  a  systematic  manner. 

The  purpose  of  this  study  is  the  application  of  modern  control 
techniques  and  optimization  procedures  to  the  design  of  a  controller 
for  a  twin  spool  variable  cycle  engine.  The  extent  of  this  study  is 
limited  to  sea  level  static  condi tiors  for  acceleration  and  decelera¬ 
tion  transients  of  differing  magnitudes.  Included  is  a  study  of  the 
controller's  ability  to  meet  steady  state  performance  and  transient- 
response  requirements  when  installation  and  deterioration  effects  are 
superimposed  on  the  nonl inear  dynamic  model . 

The  design  of  the  controller  follows  an  approach  similar  to  that 
used  by  Weinberg  (Reference  6).  In  addition,  the  original  F100-PW-100 
dynamic  model  (Reference  8)  was  modified  to  include  variable  high  and 
low  turbine  areas,  nonmixed  duct  flow  and  separate  fan  duct,  and 
core  exhaust  nozzles.  Steady  state  control  schedules  were  derived 
for  the  sea  level  static  operating  condition  similar  to  the  standard 
F100  engine.  An  interactive  computer  program  was  developed  to  solve 
the  matrix  Ricatti  equation.  The  solution  is  output  in  the  form  of 
the  feedback  gain  matrix  and  a  display  of  the  system  response  for  a 
particular  set  of  performance  index  weightings. 
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In  summary  the  objectives  of  the  study  are: 

1.  Modify  the  original  F100-PW-100  transient  model  to  a  VCE  model  and 
develop  the  steady  state  control  schedules  at  the  sea  level 
condition. 

2.  Derive  linearized  models  to  represent  the  nonlinear  model  at  the 
sea  level  condition. 

3.  Der!ve  feedback  gains  for  the  multivariable  controller  based  on 
the  linearized  models  using  an  interactive  Ricatti  equation  solver 
computer  program. 

4.  Implement  the  optimal  controller  on  the  nonlinear  model  in  a  sub- 
optimal  fashion. 

5.  Study  the  controller's  performance  with  installation  and  deterior¬ 
ation  effects. 
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SECTION  II 


VARIABLE  CYCLE  ENGINE  NONLINEAR  MODEL 
1.  VARIABLE  CYCLE  ENGINE  CONCEPT 

Reference  9  presents  a  detailed  study  of  the  benefits  of  a 
variable- turbine-geometry  turbojet  engine  versus  a  fixed-turbine-geometry 
turbofan  engine  with  comparable  capabilities.  Reference  1,  Section  3 
defines  the  engine  cycle  for  a  variable-turbine-geometry  turbofan 
engine  with  the  configuration  and  engine  station  designation  shown  in 
Figure  1.  The  increased  interest  in  variable  cycle  engines  (VCE) 
recently  is  a  result  of  requirements  for  both  military  and  super 
sonic  transport  (SST)  engines  to  provide  extremely  high  thrust  levels 
while  still  being  capable  of  operating  with  low  fuel  consumption  rates 
at  very  low  power  settings  for  subsonic  cruise. 

Currently  this  type  of  mission  profile  is  satisfied  by  an  after¬ 
burning  turbofan  engine.  The  afterburner  provides  the  required  high 
thrust  levels,  while  the  inherently  high  propulsive  efficiency  of  the  fan 
provides  good  performance  for  subsonic  cruise.  Total  aircraft  perfor¬ 
mance  suffers  for  this  engine  cycle  because  of  the  wide  variation 
of  airflows.  The  engine  cycle  requires  the  inlet  fan  airflow  to 
operate  over  a  wide  range  between  maximum  power  and  cruise  conditions. 

The  inlet  and  nozzle  must  be  sized  to  accommodate  ;.he  maximum  required 
airflow  rates,  resulting  in  high  inlet  spillage  drag  and  high  aft- 
end  drag  at  subsonic  cruise  conditions.  The  important  performance 
benefit  for  the  VCE  configuration  over  a  fixed  turbine  configura¬ 
tion  is  the  capability  to  operate  at  constant  fan  inlet  airflow 
over  not  only  the  augmented  power  range,  but  also  over  a  signifi¬ 
cant  portion  of  the  nonaugmented  high  power  range.  For  example, 
at  the  sea  level  static  condition  the  cycle  used  in  this 
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Figure  1.  Variable  Cycle  Engine  Schematic 


study  operates  at  constant  fan  inlet  airflow  from  to  100%  non- 
augmented  thrust. 

Constant  airflow  over  this  range  of  thrust  levels  is  accomplished 
through  a  mode  of  operation  referred  to  as  constant  match  varying 
temperature  (CMVT)  operation  (Reference  1).  Both  rotor  speeds  and  pres¬ 
sure  ratios  are  held  constant  as  the  high  pressure  turbine  stator  inlet 
temperature  is  varied.  As  fuel  flow  is  reduced  from  its  maximum 
value,  turbine  inlet  temperature  decreases  but  the  turbines  must 
still  maintain  a  constant  value  of  turbine  work.  This  is  accomplished 
by  varying  the  turbine  and  core  exhaust  nozzle  areas.  The  thermo¬ 
dynamic  equation  governing  this  process  is  presented  in  Appendix  A. 

This  mode  of  operation  can  continue  until  the  low  turbine  exit  flow 
parameter  reaches  a  maximum  allowable  value  as  determined  from  con¬ 
sideration  of  pressure  loss  and  flow  separation  of  the  fan  turbine 
exit  guide  vane.  Below  this  thrust  level  constant  fan  and  com¬ 
pressor  rotor  speeds  cannot  be  maintained  but  with  proper  schedul¬ 
ing  of  the  fan  duct  nozzle  and  fuel  flow  the  desired  fan  and  com¬ 
pressor  operating  lines  and  engine  bypass  ratio  can  be  controlled. 

2.  NONLINEAR  DYNAMIC  MODEL 

The  nonlinear  dynamic  model  used  for  this  study  was  derived 
from  Pratt  and  Whitney  Aircraft's  F100-PW-100  transient  engine 
simulation  (Reference  8).  This  section  will  describe  the  original 

simulation  and  the  next  section  will  describe  modifications  made  to 
represent  the  VCE  configuration  of  Figure  1.  The  engine  configura¬ 
tion  simulated  by  the  original  computer  model  is  depicted  in 
Figure  2.  The  engine  has  a  three  stage  fan  driven  by  a  two  stage  low 
pressure  turbine.  The  fan  airflow  that  bypasses  the  compressor 
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Figure  2.  Original  FlOO  Engine  Schematic  and  component 
Performance  Characteristics  (Ref.  8) 


Is  mixed  with  the  low  pressure  turbine  exit  flow  in  the  augmentor 
section  and  exhausted  to  ambient  conditions  through  a  single  variable 
geometry  balance  beam  nozzle  with  an  actuated  divergent  flap.  This 
configuration  is  known  as  a  twin-spool,  mixed  flow,  augmented  tur¬ 
bofan  engine. 

The  simulation  represents  eight  basic  components,  using  appropriate 
aerodynamic  and  thermodynamic  equations  relating  pressures,  tempera¬ 
tures,  and  mass  flows  at  various  stations  in  the  engine  and  in  terms 
of  individual  component  characteristics  as  illustrated  in  Figure  2. 

The  combination  of  the  component  characteristics,  which  are  nonlinear, 
results  in  a  nonlinear  model. 

The  calculation  flow  path  of  this  simulation  is  illustrated  in 
Figure  3-  A  steady  state  matching  point  is  always  calculated  first 
using  a  modified  Newton-Raphson  convergence  technique  (SSMITE). 
Transient  calculations  begin  at  t=0.0+At  and  iterated  until  the 
input  final  time  is  reached.  The  program  takes  about  25  seconds  Central 
Processing  (CP)  time  on  the  CDC  6600  for  a  one  second  engine  transient. 

The  engine's  dynamics  are  represented  by  the  seventeen  state 
variables  (not  including  controller  dynamics)  listed  in  Table  1. 

The  state  variables  model  the  following  three  dynamic  elements: 

1.  Torque  difference  integrations  which  set  rotor  speeds. 

2.  First-order  temperature  lags  which  represent  heat  storage. 

3.  Gas  flow  difference  integrations  which  set  pressures. 

These  states  are  controlled  by  varying  the  control  inputs;  main  fuel 
flow,  fan  and  compressor  variable  geometry,  and  exit  nozzle  area. 

The  engine  states  are  calculated  from  the  recursion  formula: 
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TABLE  T 


STATE,  INPUT,  AND  OUTPUT  VECTORS  FOR  ORIGINAL  FI 00  MODEL 
State  Vector 

Xj  =  gas  total  temperature  downstream  of  Inside  diameter  fan  exit 
guide  vanes  (TT25h) 

X2  “  total  temperature  of  gas  in  duct  (TT25c) 

x^  =  total  temperature  of  gas  at  compressor  discharge  (TT3) 

x^  =  total  temperature  of  gas  at  burner  exit  (TT4) 

xc  =  high  turbine  gas  flow  total  temperature  after  energy  loss 

!> 

to  vanes  and  blades  (TT4hi) 

Xg  =  high  turbine  gas  flow  total  temperature  after  energy  loss 
to  seals  and  discs  (TT4lo) 

Xy  -  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
vanes  and  blades  (TT**5hi) 

Xg  =  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
seals  and  discs  (TTA5I0) 

Xg  =  low  turbine  yas  flow  exit  total  temperature  (TT5) 

Xjq  =  duct  exit  gas  flow  total  temperature  (TT6c) 

Xjj  *=  augmentor  exit  gas  flow  total  temperature  (TT7) 

XI2  “  fan  physical  speed  (N1 ) 

X|g  =  compressor  physical  speed  (N2) 

x^  ■  duct  entrance  total  pressure  (PT25c) 

x^g  *=  compressor  discharge  total  pressure  (PT3) 

Xjg  88  inter-turbine  volume  total  pressure  (PT45) 

Xjy  *»  afterburner  volume  total  pressure  (PT7) 
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TABLE  1  (concluded) 


Input  Vector 

*  core  exhaust  nozzle  area  (AJC) 
u2  =  fan  Inlet  guide  vanes  (Cl VV) 
u^  =  compressor  stator  vanes  (RCVV) 
u^  =  main  burner  fuel  flow  (WF) 

u,.  **  augmentor  fuel  flow  (WFAB) 

Output  Vector 

Yj  =  net  thrust  (FN) 

y2  =  gas  temperature  at  inlet  to  high  pressure  turbine  (TTAM) 

*  airflow  through  fan  (WFAN) 
y^  ■  fan  surge  margin  (SMAF) 

y,.  »  compressor  surge  margin  (SMHC) 


(1) 


x.  .  =  X.  +t.  x  [1-e  At/,Ti,n] 
i f n+ I  i i ,n  I  ,n 

where:  x,  =  the  ith  state  after  the  nth  time  increment 
i,n 

At  -  time  increment 

=  time  constant  associated  with  ith  state  at  the  nth 


T?,n 


time  increment 


x.  =  the  time  derivative  of  the  ith  state  for  the  nth  incr'.ment 
i  >n 

The  time  rate  of  change  of  the  rotor  speeds  is  determined  by  the 
difference  between  the  shaft  torque  generated  by  the  turbine  and  the 
shaft  torque  absorbed  by  the  compressor.  These  torques  are  calculated 
from  the  turbine  and  compressor  component  thermodynamic  and  aero¬ 
dynamic  relationships.  Integration  of  the  turbine-compressor  shaft 
torque  difference  with  respect  to  time  determines  the  rotor  speed 
for  input  to  the  engine  calculations  of  Figure  3. 

Temperature  dynamics  are  caused  by  transient  heat  storage 
effects.  For  example,  a  temperature  change  at  the  inlet  to  a  turbine 
is  not  immediately  reflected  in  a  change  in  the  energy  generated  by 
the  turbine,  because  some  of  the  available  energy  goes  into  heating 
or  cooling  the  turbine  disks  and  blades. 

3.  VARIABLE  CYCLE  ENGINE  MODEL 


The  major  modifications  made  to  the  transient  simulation  described 
in  the  previous  section  included: 

1.  Changing  the  fan  duct  flow  from  mixed  flow  to  nonmixed  flow. 
This  necessitated  the  addition  of  a  fan  duct  exhaust  nozzle. 

2.  High  and  low  pressure  turbine  flow  parameter  and  efficiency 
multipliers  to  simulate  the  effect  of  variable  turbine  geometry. 

3.  Elimination  of  afterburner  logic. 


12 


4.  Elimination  of  control  scheduling  logic  and  actuator 
dynamics. 

5.  Addition  of  steady  state  control  schedules  for  the  sea 
level  flight  condition  used  for  this  study. 

The  VCE  concept  has  previously  been  discussed  and  will  not  be 
repeated  in  this  section.  To  change  the  fan  duct  to  a  nonmixed 
flow  configuration  required  a  separate  exhaust  nozzle  for  the  flow 
that  bypasses  the  compressor.  The  fan  duct  nozzle  calculations  were 
implemented  similar  to  those  for  the  core  nozzle.  Both  the  core 
and  fan  duct  exhaust  nozzles  are  variable  and  are  represented  by 
discharge  coefficients  and  velocity  coefficients.  The  static  pressure 
balance  at  Station  6  of  Figure  1  required  by  the  original  model  to 
have  a  balanced  steady  state  condition  was  eliminated.  The  elimination 
of  this  allows  the  fan  operating  point  to  be  constant  in  the  CMVT 
range  of  operation. 

The  addition  of  variable  area  turbines  was  accommodated  through 
multipliers  on  efficiency  and  flow  parameter  as  a  function  of  vane 
angle.  These mul tipi  iers  are  presented  in  Figure  4.  It  was  assumed 
that  the  turbine  design  point  was  at  intermediate  power  (PLA  =  83). 
There  was  no  unilateral  turbine  efficiency  loss  assumed  from  the 
original  model.  This  would  most  likely  not  be  true  for  an  actual 
case  but  the  effect  on  this  study  was  considered  minimal.  The 
relationships  of  Figure  4  are  simplified  best  estimates  from  the 
information  available.  No  actual  data  was  available  for  the  cycle 
configuration  used  in  this  study. 

The  scope  of  this  study  did  not  include  the  use  of  after¬ 
burning  which  was  simulated  in  the  original  model.  Therefore,  the 
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afterburner  logic  was  eliminated  to  reduce  computer  memory  requirements. 
The  original  control  scheduling  logic  and  actuator  dynamics  were  also 
eliminated  to  reduce  the  memory  size.  The  VCE  simulation  derived  for 
this  study  did  not  include  actuator  dynamics.  This  study  represented 
a  preliminary  design  effort  for  which  che  act'ator  and  sensor  responses 
were  assumed  to  be  instantaneous.  The  control  designed  in  this  study 
did  not  move  the  actuators  instantaneously. 

The  steady  state  control  schedules  derived  for  this  study  are  shown 
in  Figure  5.  These  schedules  are  applicable  only  at  sea  level  static 
and  were  derived  based  on  certain  constraints.  For  the  CMVT  power  range 
the  fan  and  compressor  rotor  speeds  and  pressure  ratios  were  held  to  the 
same  values  of  the  original  model  operating  at  intermediate  power  (PLA  = 
83).  This  resulted  in  the  fan  duct  exhaust  nozzle  area  (AJD)  and  fan 
and  compressor  vane  position  (CIVV,  RCW)  being  constant  in  the  CMVT 
power  range.  The  relationship  between  fuel  flow  (WF)  and  power  level 
angle  (PLA)  was  arbitrarily  established  to  be  approximately  the  same  as 
in  the  original  model.  The  two  turbines  and  core  nozzle  area  schedules 
in  the  CMVT  mode  were  _’-*.termined  using  SSMITE.  The  three  variable  areas 
were  allowed  to  vary  until  the  fan  and  compressor  were  operating  at  the 
proper  mstch  point.  These  areas  were  calculated  for  fuel  flows  from 
intermediate  power  to  the  CMVT  breakpoint.  The  CMVT  breakpoint  occurred 
at  PLA  =  55,  as  determined  from  the  considerations  mentioned  previously. 

For  power  settings  below  the  CMVT  breakpoint  the  two  turbines  and 
core  nozzle  areas  were  held  constant  and  fan  and  compressor  rotor 
speeds  fell  off.  The  prescribed  fan  operating  line  was  maintained 
by  allowing  the  fan  duct  nozzle  to  vary  and  power  was  reduced  by 
reducing  fuel  flow.  The  fan  and  compressor  guide  vanes  were  initially 

15 


scheduled  as  a  function  of  their  respective  rotor  speed  using  the 
schedules  from  the  original  model.  Once  the  speeds  were  determined  as 
functions  of  PLA,  the  CIVV  and  RCVV  schedules  were  made  a  function  of 
PLA  instead  of  speed.  This  change  in  the  scheduling  function  was 
necessitated  by  the  methodology  used  to  derive  the  multivariable 
control  law  (this  will  be  explained  in  a  later  section).  From  a  steady 
state  performance  standpoint  it  would  be  better  to  schedule  the  guide 
vanes  as  a  function  of  rotor  speed. 

These  steady  state  schedules  make  it  possible  to  run  the  simu¬ 
lation  at  any  power  setting  (PLA)  between  20  and  83  and  obtain  steady 
state  output.  Figure  6  presents  plots  of  some  of  the  output  and  state 
variables  vs.  %  intermediate  thrust.  The  rest  of  tnis  report  will 
discuss  a  multivariable  control  that  controls  the  simulation  during 
large  or  small  transients  over  the  entire  PLA  range. 


ermediate  FN 


SECTION  III 


CONTROL  DESIGN  METHODOLOGY 
1.  LINEARIZATION  TECHNIQUE 

The  variable  cycle  engine  used  in  this  study  is  represented  by 
the  18  state  variables  presented  in  Table  2.  This  section  summarizes 
a  linearization  technique  developed  by  Weinberg  (Reference  10)  which 
linearizes  a  high  order,  nonlinear  digital  simulation  to  state  vari¬ 
able  format  by  taking  partial  derivatives  of  the  state  variable  time 
derivatives  with  respect  to  the  states  and  inputs.  The  system's 
eigenvalues  are  then  investigated,  as  described  in  the  next  section, 
and  a  reduction  is  performed  so  that  only  the  effects  of  dominant 
poles  remain. 

Optimal  control  theory  is  based  on  a  mathematical  model  in  state 
variable  form.  The  state  variable  form  is  convenient  because  the 
differential  equations  are  ideally  suited  for  digital  solution,  it 
provides  a  unified  framework  for  the  study  of  nonlinear  and  linear 
systems  and  it  has  strong  physical  motivation  (Reference  11).  The  lin¬ 
earization  converts  the  eighteen  first  order  nonlinear  equations  of  the 
VCE  simulation  indicated  by  Equation  1  to  a  linear,  time  invariant 
state  variable  format  (Reference  10); 

x  =  Ax  +  Bu 

y  =  Cx  +  Du 

where;  x  -  eighteenth  order  state  vector 
u  =  seventh  order  input  vector 
y  «  fifth  order  output  vector 
The  overline  indicates  a  vector  or  matrix,  the  overdot  represents  a 
time  derivative.  A  is  the  (18  x  18)  matrix  of  partial  derivatives  of 


(2) 

(3) 


TABLE  2 


STATE,  INPUT,  AND  OUTPUT  VECTORS  FOR  VCE  MODEL 
State  Vector 

Xj  a  gas  total  temperature  downstream  of  Inside  diameter  fan  exit 
guide  vanes  (TT25h) 

x^  -  total  temperature  of  gas  in  duct  (TT25c) 

x^  =  total  temperature  of  gas  at  compressor  discharge  (TT3) 

x ^  =  total  temperature  of  gas  at  burner  exit  (TT4) 

xc  =  high  turbine  gas  flow  total  temperature  after  energy  loss 

b 

to  vanes  and  blades  (TT4hi) 

Xg  =  high  turbine  gas  flow  total  temperature  after  energy  loss  to 
to  seals  and  discs  (TT4lo) 

x^  =  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
vanes  and  blades  (TT45hi) 

Xg  “  low  turbine  gas  flow  total  temperature  after  energy  loss  to 
seals  and  discs  (TT451o) 

xG  =  low  turbine  gas  flow  exit  total  temperature  (TT5) 

Xjq  =  duct  exit  gas  flow  total  temperature  (TT6c) 

Xjj  -  augmentor  exit  gas  flow  total  temperature  (TT7) 

Xj2  a  fan  physical  speed  (Nl) 

Xj^  "  compressor  physical  speed  (N2) 
x^  “  duct  entrance  total  pressure  (PT25c) 
xlc  =  compressor  discharge  total  pressure  (PT3) 

Xjg  85  inter-turbine  volume  total  pressure  (PT45) 

Xjy  **  afterburner  volume  total  pressure  (PT7) 
x^g  «  fan  duct  volume  total  pressure  (PT6c) 
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TABLE  2  (concluded) 
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1C 


the  state  time  derivatives  with  respect  to  the  states; 
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C  is  the  (5x1 8)  matrix  of  partial  derivatives  of  the  outputs  with 
respect  to  the  states; 


3Y 

3x 


1 

1 


ax, 8 


c  = 


•  •  •  •  • 

3x,  3J£,8 


(6) 


5  is  the  (5x7)  matrix  of  partial  derivatives  of  the  outputs  with 
respect  to  the  inputs; 


3y, . 3y?  | 

^  ^  i 


3y5 . 3y5 

3u ,  3Uj 


The  basic  strategy  used  to  calculate  the  state  matrices  of 
Equations  4  thru  7  was  to  converge  on  a  steady-state  condition  and 
then  perturb,  one  at  a  time,  each  of  the  state  and  input  variables. 
The  nonlinear  dynamic  simulation  was  modified  to  incorporate  the 
following  algorithm  for  calculating  the  linearized  state  matrices; 
1.  A  steady  state  point  at  PLA  =  83  is  run  to  obtain  the 
normalizing  values. 


2.  A  steady  state  point  at  the  desired  power  setting  is  run  to 


obtain  x,  (the  choice  of  power  settings  will  be  discussed 
i » n 

later) . 

3.  Each  of  the  state  and  input  variables  (25  in  all)  are 

perturbed  one  at  a  time  from  their  steady-state  values  and 
the  transient  simulation  run  for  fifty  time  steps  to 
obtain  x.  n+j  of  the  A  and  B  matrices  in  Equations  k  and  5. 

k.  The  C  and  D  matrices  of  Equations  6  and  7  are  calculated  in 
a  similar  manner  as  for  A  and  B;  however,  no  integrations  were 
required  to  perform  the  derivative  calculations  of  Equations 
6  and  7* 

The  derivatives  were  calculated  from  50  time  steps  of  0.00001 
seconds  each.  The  state  and  input  variables  were  usually  incremented 
0.2%  of  their  base  values.  Once  calculated  these  matrices  were 
punched  out  on  cards  in  a  format  suitable  for  input  to  the  matrix 
reduction  program  discussed  in  the  next  section.  For  more  specific 
details  of  the  above  procedure  see  Reference  10. 

2.  MATRIX  REDUCTION 

To  facilitate  the  analysis  and  design  of  dynamic  control  systems, 
it  is  often  desirable  to  reduce  the  system's  order.  In  linear 
systems,  this  reduction  is  possible  when  the  system's  transfer  func¬ 
tion  contains  dominant  poles  which  determine  the  transient  response 
essentially  independently  of  the  remaining  poles  (Reference  12).  In  the 
complex  plane,  the  less  dominant  poles  appear  to  the  left  of  the 
dominant  poles  and  represent  transient  terms  which  die  quickl>  when 
compared  to  the  transient  terms  associated  with  the  dominant  poles. 

This  section  summarizes  a  technique  developed  by  Weinberg  (Reference  12). 
which  reduces  a  high  order  system  to  a  lower  order  system,  whose 


poles  correspond  to  the  dominant  poles  of  the  higher  order  transfer 
function.  This  reduction  technique  does  not  require  previous 
knowledge  of  system  characteristics  and  identifies  important  dynamic 
states  in  a  straight-forward  manner. 

The  values  of  the  A,  B,  C  and  5  for  the  1 8th  order  linearized 
modes  are  given  in  Table  3.  An  algorithm  (References  10  and  12)  to 
reduce  the  order  of  the  model  is  summarized  below. 

Recalling  the  state  variable  form  of  Equations  2  and  3  it  is 
assumed  that  the  matrix  A  has  r  nondominant  poles  which  can  be 
eliminated,  implying  that  r  states  of  x  have  zero  derivatives. 


x  =  Ax  + 

Bu 

(8) 

y  =  Cx  + 

Du 

(9) 

Equations  8  and  9  can  be  rewritten  with  partitioned  matrices; 


Jel 

Apll  Apl2 

+ 

B  . 

-Bi. 

• 

V 

A  A 

p21  p22 

Xp2 

B  9 

P2 

(10) 


y  = 


Cp1  Cp2 


.il. 

XP2 


+  D  u 


(11) 


where:  x  .  =  (n-r)xl  vector  whose  derivatives  x  ,  are  dynamically 
pi  pi 

important 

•  £ 

Xp^  *=  (rxl)  vector  whose  derivatives  x^2  are  dynamically 

i nsignif leant 

(n-r)xl 

(n-r)xr  matrix 


A,,3  (n-r)x(n-r)  matrix 
pi  1 


pi  2 

Ap2i  3  rx(n-r)  matrix 
»  (rxr)  matrix 


TABLE  3 


18th  Order  Linearized  Matrices  for  PLA=60 
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23.366Ju 

3o  6 «  0  b  y  3  0 

-546.  *>0  7  C  0 

VARIA -ILL  17 

-.06770 

D. CuOCC 

-.  o.’SaO 

-2. J7570 

•5o260 

. 06  3J  3 

-.00020 

-  0. 3000  0 

.03370 
-7  •  lo  190 

rt  «  0  C  0  C  j 
•  UC  li  *•  c 

2.  ‘.-19  '  J 

-.23654 

•  1-.910 

*  337>0 

•  76 150 

*  3.79-60 

VARIA  9LE  13 

12.10120 

. 112  aO 

-4  0?iU v 

•  02800 

-.00790 

-.00090 

.00750 

.  j  J  C  8  0 

-0.03330 
-.00  9-0 

3.  JJl’O 
-:*.  71330 

-  0  *  COOOi 

12.03220 

-13. 77-00 

-7 .03930 

• C 151 0 

.0194.0 
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TABLE  3  (continued) 


P  MATRIX 
VARlA-Kf  1 

r.ocooo 

-.m.iu 

.  Ol-’T 

-0. 31030 

VAPIA  >26  2 

-e.  coioo 

-  •  13-  jU 

-. CO" K 

0. 03030 

VtKIAPU  1 
-o.ooooo 

.00  370 

-  •  2  !  6 1 2 

.00-120 

VASIA  SUI  1. 
-c.oonOb 

.31b"0 

A.6--0C 

23. OlaR G 

VA«1A3l<T  3 
C.3J009 

.00330 

.  07O-33 

.*-2  O-.C 

VARlfOLS  o 

C. 00OJ0 

C.  JOliO 

. 0011- 

.  0)7’0 

VARifur  r 
o.Ougo: 

.0C-.20 

. 06322 

lillsU- 

VAkIAUI.  5 

n.  00:00 

0 .3C-3  '33 

. 03272 

.  9-  620 

VASIIEL5  9 
-.  00  320 

. 00320 

.b>7  31. 

.66310 

VA-ilMLC  10 
-C.CC20  0 

.0331 3 

-.117-. 

.  OObCC 

VAi-'IA  OLE  11 

-c.cj  no 

0.03303 

.u3-’: 

. G-2i0 

VARIA3LK  12 
•  OU  170 

.32?'  3 

-» Cl°-3 

-.01603 

VAR  It  OLt*  13 
-0.00330 

.oer  c 

.  boJ  30 

- .  20 '  3C 

VARKeLC  H. 
O.OOOJ3 

-.so  ja: 

1. 721‘,1 

-0.3jo-’C 

VARIASlF  Is 
0.8G099 

-.13.--0 

oO  o’  * 

3.  15393 

VAftXAJlE  16 

0.0020b 

-  .  G22C0 

-13.  U633 

5.157’f 

VASIA -3LF  17 
-2. -1593 

-.ocau 

-.ci.-n: 

-.  02750 

VAR1A  OIL  10 
C.C0  300 

-1.767- 0 

1.  G716C 

.  3.1C20 

C  MAT  P  IX 
VA-iA'IlF.  l 
-.00  13  0 

0.00333 

C. C300C 

0.  53  3  3C 

C. 30 00 3 

-.Cb2  -.0 

-.03*.  13, 

0.  03C3C 

1 . 0-4  33  3 

VAPJAUE  2 
•  co  *2  2 

.392 '"3 

1) .  3  C  3  •  3 

.Otrl- 

.97210 

0.00333 

0.30-10 

0.  C 3 r  3. 

-0. 30202 

0 • n0  330 
VAR1A3LC  3 
-.30013 

.3. 00  3  '  -3 

3. 9t  300 

O.OOOO: 

C. 00330 

C.  03  303 

■-.  300  3 

J.  000 17 

.?99oC 

0.00030 
VARIA  ?lt  - 
-.01130 

0.  )C3-1 

0.3C-30  j 

C. COO"0 

0.  OC  090 

0.30-3  33 

o.ocrco 

0.  U33  33 

1.7 lu5G 

0  .  Ca  JOb 
VAOIA'iU.  6 
-.09  270 

o.ouoro 

r.ico'o 

o.ic-::. 

0.0)200 

0.00000 

0*0033: 

e.  30Cj: 

.90s  30 

r. 00330 

P  MATRIX 

VAR  11  •'’1C  1 

.29370 

0.  iOOCO 

-0.30  ot  3 

-.00020 

-. 30130 

VARIA-IU'.  2 

2. H003C 

P  •  u  0  0  0  0 

.  C096-: 

• UC-lbb 

VARIA  lit  3 

o.cooo: 

-.1C.3.-0 

-(1.0003. 

0. 03203 

VAKiAELC  A 

C.  00-300 

-.12330 

102-J.I 

C. 03330 

VARiA'Rf  3 
0.00300 

-.31.0  !o 

-.2ji  33 

u «  00  u  ;  0 

—  3. 00033 

-.30010 

-.000-0 

- .  30013 

-  3 . 0  3  3  0  3 

.00310 

0.0J0J0 

.  u  3  7  6  0 

-. OOClO 

.0u020 

.00-30 

. b022C 

3.30330 

•  .0013 

0.00002 

0  *  20000 

0.00000 

0. C00J0 

-3.03300 

-.JOs-iO 

.29030 

-3.03000 

-.  303'*3 

.01170 

0*00033 

.0696) 

-.170-0 

-  3 .0 3u30 

-.50310 

O.bOOlO 

0.23303 

.00173 

-.00030 

0.03000 

-  .  1 . 3  0  0 

-.27-20 

0.00033 

-.7-5 lj 

-.b-960 

-  •  196  7  u 

.33 1-2 

-.00033 

-.131-0 

3  0 . 7  0  1  f.  3 

-. 001-u 

-.900-0 

-llo.bo2D3 

113.3—  70 

-0.32030 

- .  1--30 

- • 7  o73  0 

-6.60930 

.33250 

- .  0  2  b  3  3 

0.30233 

,0.u*030 

3 . 0  3  0  3  0 

0  •  j  3  3  0  0 

.03113 

« 0  0  OoO 

-u .000)9 

0 .33:20 

0  *  3- u 0 0 

J  .  0  j  0  )  0 

3. 00030 

3. 3  0  *,  2  0 

.0  5r-l0 

-.00130 

3 . b  0  G  3  0 

0.33350 

3.33003 

o.:j  }  ] 0 

3 . 0  0  b3C 

0.  JC  JCO 

.0 j020 

-  •  0  0  ?6t) 

-3.  0003b 

0, 30  30  0 

0. 30300 

0  •  0  0  C  0 1> 

3. 00 030 

C. 9019  0 

.01230 

-  •  «>*•  32  2 

-0. 00030 

C.  2 0  3c  0 

3.03000 

O.uOCGO 

O.OOolO 

0. 9  3  J  0  0 

i  • T  0  2  2  3 

•  ^2  J 

-, 797  J  0 

0.33300 

.20230 

0  •  J  0  0  9  0 

O.OObDO 

0.30000 

0.90000 

0.  0  0  33C 

3.00009 

0.30093 

9, 0003C 

0,00030 

3.000)3 

0.00030 

0.0003b 

0.20003 

■3.  2993b 

B  ,  »  (n-r)xm  matrix 
P« 

B  0  -  (r-m)  matrix 
P2 

Cpj  ®  px(n-r)  matrix 

Cp2  *•  (pxr)  matrix 

D  =  (pxm)  matrix 

n  =  number  rf  states  (1 8) 

m  =  number  of  control  inputs  (7) 

p  =  number  of  outputs  (5) 

r  =  number  of  non-dominant  poles 

Reference  j 2,  Appendix  A  presents  an  algorithm  for  converting  a  set  of 
state  equations  into  the  form  of  Equations  10  and  II.  Setting  xp2  =  0, 
the  reduced  state  equations  become; 


*p1  =  Vpl  +  §r“  (,2> 

y  "  Crxpl  +  Dru  (13) 

where:  Af  -  ApH  -  Apl2Ap22  lAp21 

^r  =  ®p1  "  \u\ll  1gp2 

^r  =  Ep1  "  ^p2%22  ^p21 

^r  =  ^  "  Cp2Ap22  1 ®p2 

Equations  12  and  13  completely  define  the  reduced  system  model 
but  the  state  derivatives  of  x  _  must  be  determined  before  the  reduc- 

p2 

tion  can  be  performed.  The  eigenvalues  for  the  VCE  system  are  given 
in  Table  4.  Eigenvalues  with  real  parts  less  than  -15  are  considered 
to  be  nondominant  poles.  It  is  observed  that  the  system  derived  in 
this  study  has  six  dynamically  important  poles. 

Summarizing,  to  determine  the  unimportant  state  derivatives  of 
Xp2  the  following  procedure  was  used; 
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TABLE  4 

EIGENVALUES  OF  A  MATRIX  FOR  PLA  =  60 


18th  order 

6th  order 

4th  order 

-547.08 

-12.27 

-12.34 

-130.03 

”0.66 

-3. 66+1. 57 j 

-64.68 

-3.70+1 .55j 

-3.65 

-57.12 

-3.64 

-50.92 

-2.00 

-49.38 

-36.24 

-24.47i-6.45j 

-22.54 

-18. 58+3 • 53j 

-9.56 

-3.65il.8lj 

-3.04 

”0.66 

-1.99 


1.  Assume  that  a  pole  (or  poles)  Is  nondominant. 

2.  Rearrange  x  into  the  submatriceS  Xpj  and  x^* 

3.  Rearrange  A,  B,  and  C  to  correspond  to  the  new  order  of  x. 

4.  Calculate  the  reduced  matrix  Ar  and  determine  its  eigenvalues. 

5.  If  x^  consists  of  only  derivatives  corresponding  to  non- 
dominate  poles,  the  dominant  poles  of  Ar  will  approximate 
the  dominant  poles  of  A. 

A  program  from  Referenced  was  used  to  automatically  perform  steps  2 

thru  4.  The  calculated  eigenvalues  for  a  sixth  order  system  ard  a  fourth 

order  system  are  presented  in  Table  4.  The  states  used  in  the  submatrix 

Xpj  for  the  sixth  and  fourth  order  systems  are  presented  in  Tible  5. 

These  are  the  states  that,  when  individually  substituted  as  x  „  into 

p2 

the  matrix  reduction  program,  caused  the  dominant  poles  to  shift.  The 
resulting  A,  B,  C,and  6  matrices  for  the  sixth  and  fourth  order  systems 
are  presented  in  Table  6. 

The  reduction  from  six  to  four  states  was  possible  because  two  of 
the  states  had  no  effect  on  the  output  variables  as  evidenced  by  the 
two  columns  of  near  zero  values  in  the  C  matrix,  6th  order,  of  Table  6. 
This  is  further  demonstrated  in  Figure  7  where  the  output  variables 
of  the  1 8th,  6th,  and  4th  order  linear  models  and  the  output  of  the 
nonlinear  model  are  plotted  for  a  step  Input  in  fuel  flow.  The  actual 
engine  is,  in  this  case,  modeled  W  18  states  that  denote  energy 
storage.  The  four  states  identified  by  the  reduction  program  are  the 
states  which  store  energy  sufficient  enough  to  affect  the  engine's 
transient  performance.  The  fourth  order  model  at  two  different  power 
settings  was  used  for  this  study  to  derive  a  multivariable  control  as 


described  in  inter  sections. 


TABLE  5 


STATE  VECTORS  FOR  6TH  AND  4TH  ORDER  MODELS 

6th  order 

Xj  TT^llo  -(x^  )-  high  turbine  gas  flow  total  temperature  after  energy 

loss  to  vanes  and  blades 

x2  TT451o  -(x8  )  -  low  turbine  gas  flow  total  temperature  after  energy 

loss  to  seals  and  discs 


x3 

m 

~  (x 1 2 )  ~  ^an  phy$5cal  sPeec^ 

x 4 

N2 

-(x^)_  compressor  physical  speed 

*5 

PT7 

-(xj^)-  afterburner  volume  total  pressure 

x6 

PT6c 

-(xlg)-  fan  duct  volume  total  pressure 

4th  order 

X1 

N1 

-(xj2)-  fan  physical  speed 

x2 

N2 

compressor  physical  speed 

x3 

PT7 

- (x j y) -  afterburner  volume  total  pressure 

x4 

PT6c 

-(xjg)-  ^an  duct  vo^ume  total  pressure 

(  )-  state  variable  number  from  Table  ; 

Vt.fi  It  «LF 


3.  CONTROL  OPTIMIZATION  TECHNIQUE  FOR  LINEARIZED  MODELS 


As  stated  in  the  introduction,  the  design  of  a  multi-input, 
multi-output  control  system  using  classical  techniques  requires  a 
significant  amount  of  trial -and-error,  experience  and  intuition.  The 
purpose  of  this  study  was  to  further  demonstrate  the  application  of 
optimization  techniques  as  a  design  tool  for  multivariable  control 
systems.  This  section  will  discuss  some  of  these  optimization  tech¬ 
niques  and  their  application  to  the  control  design  for  the  VCE. 

The  objective  of  optimal  control  theory  is  to  determine  the 
control  signals  that  will  cause  a  process  to  satisfy  the  physical 
constraints  and  at  the  same  time  minimize  some  performance  criterion 
(Reference  11).  For  the  variable  area  turbine  engine,  some  of  the  phys¬ 
ical  constraints  include  fan  and  compressor  surge  marg ins, maxi  mum  turbine 
inlet  temperature  level  and  rate  of  change,  maximum  pressure  levels, 
maximum  and  minimum  fuel-air  ratios  for  burner  stability  and  the  rate 
and  amplitude  limits  of  all  control  variables  (Reference  7).  These  con¬ 
straints  must  be  observed  both  in  steady  state  and  in  transient  modes. 

In  addition  to  these  constraints  there  are  certain  performance  require¬ 
ments  such  as  providing  correct  steady  state  thrust  and  airflow  levels 
and  transient  thrust  response  times.  This  set  of  constraints  and 
requirements  is  independent  of  the  method  used  to  design  the  control. 
Although  implementation  as  a  controller  may  be  impractical,  the 
optimal  control  law  provides  insight  useful  in  designing  a  more  easily 
implementable  suboptimal  controller. 

A  special  class  of  optimal  control  theory  is  the  linear  quadratic 
regulator  (LQR)  problem.  For  the  LQR  problem  applied  to  a  turbine 
engine,  the  optimal  control  law  is  a  linear  time-invariant  function  of 


the  system  states.  For  the  general  system  described  by  the  linear 
state  equations;  (Reference  *  . 

x(t)  =  A(t)x(t)  +  B(t)u(t)  0*0 

The  performance  index  (PI)  to  be  minimized  is; 

PI  =  h  xT(tf)Hx(tf)  +hfJ  [xT(t)Q(t)x(t)  +  u(t)R(t)u(t)]dt  (15) 
r  r  c0 

where:  t^  =  final  time  and  is  fixed 

H  and  Q  are  real  symmetric  positive  semi-definite  matrices 
R  is  real  symmetric  positive  definite  matrix 
It  is  assumed  that  the  states  and  controls  are  not  bounded.  The 
physical  significance  attached  to  the  performance  index  is:  It  is 
desired  to  maintain  the  state  vector,  x,  close  to  the  origin  with¬ 
out  an  excessive  expenditure  of  control,  u,  effort. 

For  the  situation  wherein  the  process  is  to  be  controlled  for 
an  infinite  duration,  t^->-  it  can  be  shown  (Reference  11)  that  if  the 

plant  is  completely  controllable  and  II  =  0  than  A,  B,  R,  and  Q  are 
constant  matrices.  Equations  \k  and  15  then  become; 

x(t)  =  Ax(t)  +  Bu(t)  (16) 

PI  =  V*  [xT(t)Qx(t)  f  uT(t)Ru(t)]dt  (17) 

0 

The  optimal  control  law  for  this  case  is  derived  using  the 
Hamiltonian  ^Reference  11); 

u*(t)  =  -R*lBTKx(t)  (18) 

or 

G*(t)  =  mo 

where  K  is  obtained  by  solving  the  nonlinear  Ricatti  equation; 

0  =  -KA  -  ATK  -  Q  +  KBR'1^  (19) 

Before  proceeding  with  the  solution  of  the  matrix  Ricatti  equation 


It  Is  necessary  to  rewrite  Equations  16  and  17.  Recalling  the  discus¬ 
sion  on  the  linearization  procedure,  x  and  u  were  actually  normalized 
deltas.  For  example  if  x^  was  fan  speed,  N I ,  then  the  linearization 

would  have  used  ANI  =  (Nl  -  Nl  )/Nl  .  Where  N1  is  fan  speed  at 

ss  p  n.  ss  r 

steady  state,  Nl^  is  perturbed  fan  speed  from  steady  state  and  Nl^  is 
the  fan  speed  at  PLA  =  83  used  to  normalize.  Therefore^  Equation  16 
should  be; 

(Ax(t))  =  A(Ax(t) )  +  §(Au  (t) ) 


To  simplify  the  derivations  that  follow,  the  A's,  t's  and  overlines 
will  be  dropped. 

Next,  recall  that  the  linear  system  is  described  by: 


x  =  Ax  +  Bu 
y  =  Cx  +  Du 

let  X  =  1  xl  and  U  =  u 

I  u  I 


• 

A  1 

B 

Q 

then  X  = 

n 

0 

X  + 

1 

(20) 

(21) 


or 


X  = 

X 

< 

+  B 

'U 

Y  = 

y 

SS 

y 

X 

u 

h 

T2 

n 

^2  “  2  ^2  i 

1=1 


n 

[  E 


(22) 


Where  </^  is  the  rate  of  change  of 
Table  2  and  calculated  according  to  the 


following  formula; 
m 


m 


B.  .u.]  +  E 


>1 


A.  .x.  +  E 
i.J  J  j=1  »tJ  J 


i=1 


D0  .U. 

2,i  i 


C  1  D 

_0_ 

Y  = 

1 ! 0 

X  + 

0 

0  1  1 

"5" 

C^jc'B 

■ 

1  CM 

1  Q 

Y  = 

CX  +  D 

u 

(23) 


(2*0 
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The  parameters  of  X  and  Y  are  presented  In  Table  7  for  the  VCE 
nodel.  The  redefined  state  vector  now  Includes  the  original  states 
and  the  control  input  variables.  The  redefined  output  vector  includes 
the  original  output  variables  plus  the  redefined  state  vector  and 
the  time  derivative  of  turbine  inlet  temperature.  The  control  vector 
was  redefined  to  be  the  rate  of  change  for  the  control  inputs-  A 
FORTRAN  program  that  automatically  transforms  the  original  A,  B,  C  and 
D  matrices  into  the  form  of  Equations  22  and  2k  was  developed. 

The  purpose  of  these  redefinitions  was  to  allow  the  performance 
index.  Equation  17,  to  be  rewritten  to  allow  the  designer  direct 
weighting  of  the  state,  output,  control,  and  control  rate  variables. 
This  is  accomplished  by  adjusting  the  weighting  matrices,  Q  and  R 
in  Equation  17.  Adjusting  the  relative  weights  alters  the  solution 
of  the  minimization  such  that  the  more  heavily  weighted  terms  are 
driven  to  zero  more  rapidly  than  the  lightly  weighted  terms. 


The  desired  form  for  the  performance  inde '<  is; 


PI  »  Hf  [YTQY  +  UTRU]dt 

0 

(25) 

and  the  optimal  control  law  is; 

*  -IT 

U  =  -R  B*KX 

(26) 

This  is  no  longer  in  the  LQR  format.  Equation  25  is  transformed  by 
substituting  for  Y  with  Equation  2k\ 

YTQY  =  [CX  +  D'U]TQ[C'X  +  D'U] 
this  is  equivalent  to; 

YTQY  »  [  C'X]T(i[C'X]  +  [D'U]TQ[D'U]  +  [C'X]TQ[D'U]  +  [D'U]TQ[C'X] 
which  reduces  to; 

ytqy  =  xtc,tqc'x  +  utd,tqd'u  +  XTC,TQD'U  +  UTD,TQC'X 


REDEFINED  STATE,  INPUT,  AND  OUTPUT  VECTORS 


State  Vector 

Xj  =  fan  physical  speed  (N1 ) 

X2  =  compressor  physical  speed  (N2) 

X^  »  afterburner  volume  total  pressure  (PT7) 

X^  *»  fan  duct  volume  total  pressure  (PT6c) 

Xc  =  core  exhaust  nozzle  area  (AJC) 
i> 

Xg  =>  fan  inlet  guide  vane  angle  (CIW) 

Xj  **  compressor  stator  vane  angle  (ROW) 

Xg  =  main  burner  fuel  flow  (WF) 

X-  «  fan  duct  exhaust  nozzle  area  (AJD) 

9 

XJQ  »  high  turbine  nozzle  position  (UTVPOS) 
X^  **  fan  turbine  nozzle  position  (FTVPOS) 


Input  Vector 


U1 

U2 

U3 

Ui» 

U5 

U6 

U-, 


time  derivative 

II  II  II 

II  II  II 

II  II  II 

It  II  II 

II  II  II 

II  II  II 


of  core  exhaust  nozzle  area  (AJC) 
fan  inlet  guide  vane  angle  (CIW) 
compressor  stator  vane  angle  (RCVV) 
main  burner  fuel  flow  (WF) 
fan  duct  exhaust  nozzle  area  (AJD) 
high  turbine  nozzle  position  (HTVPOS) 
fan  turbine  nozzle  position  (FTVPOS) 


TABLE  7  (concluded) 


Output  Vector 

Yj  =  net  thrust  (FN) 

Y2  =  gas  temperature  at  inlet  to  high  pressure  turbine  (TT^M) 

Y^  =  airflow  through  fan  LWFAN) 

Y^  =  fan  surge  margin  (SMAF) 

Y,.  =  compressor  surge  margin  (SMHC) 

Y^  ■  fan  physical  speed  (N1) 

Y^  “  compressor  physical  speed  (N2) 

Yg  =  afterburner  volume  total  pressure  (PT7) 

Ya  85  fan  duct  volume  total  pressure  (PT6c) 

3 

*=  core  exhaust  nozzle  area  (AJC) 

Yj  j  13  fan  inlet  guide  vane  angle  (IGV) 

Y^2  =  compressor  stator  vane  angle  (RCW) 

Y^  **  main  burner  fuel  flow  (WF) 

Yj^  3  fan  duct  exhaust  nozzle  area  (AJD) 

Y,  =  high  turbine  nozzle  position  (HTVPOS) 

Yjg  -  fan  turbine  nozzle  position  (FTVPOS) 

Yjy  -  time  derivative  of  gas  total  temperature  at  inlet  to 
high  turbine  (TTAM) 
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and  Equation  25  becomes; 
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According  to  Reference  \k  (pg.157)  Equation  27  is  equivalent  to; 
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dt 
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By  letting;  R'  »  Q^+R  and  Q'  -  *Q2  Eqi  '/.6  (...‘duces  to 

°°  T  T 

PI  =  =5/  [X'Q'X  +  U  R'U]dt  (29) 

0 

The  optimal  control  law  is  given  by; 

U*  «  R,H{Q2  +  B'TK)X  (30) 

or 

U*  -  F'X  where;  F'  =  R'”1 (qJ  +  B'TK) 

Where  K  is  obtained  by  solving; 

0  =  — K(A*  — BR* "”1  qJ)  "  (A* -B'R,_1qJ) TK  +  KB'R,“1B,TK  -  (Qj-^R'^qJ) 
which  reduces  to; 

0  ■  -KA'  -  A'TK  +  KB'R,-1Q2  -  a' 

in  summary  the  important  equations  for  the  quadratic  requlator 
problem  as  used  in  this  report  are; 
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4.  INTERACTIVE  RICATTI  EQUATION  SOLVER 

A  FORTRAN  program  that  solves  the  Ricattl  equation  and  outputs 
the  feedback  gain  matrix  of  Equation  30  was  developed  from  a  program 
used  in  Reference  6.  The  linear  matrix  subroutines  are  from  Reference  13. 
The  program  was  written  for  a  CDC  6600  computer  to  be  operated  inter¬ 
actively  using  a  Tektronix  4014-1  interactive  graphics  unit.  The  user 
must  supply  a  set  of  A,  B,  C,and  D  matrices  from  Equations  22  and  24, 
a  set  of  state  initial  conditions,  X,  and  a  set  of  values  for  the 
diagonal  matrices  Q  and  R.  With  this  information  the  program  solves 
for  K  in  Equation  31  and  calculates  the  feedback  gain  matrix  which  is 
displayed  and  optionally  punched  out  for  use  in  a  controller  on  the 
nonlinear  transient  simulation.  Additionally  the  proqram  displays 
on-line,  transient  responses  of  all  the  output  v'r  ab'c.'  in  Table  7- 
The  optimal  control  law  of  Equation  30  is  opt.*M1  <  nl>  for  the 
particular  set  of  Q  and  R  matrices  input  to  the  program  and  within 
the  linear  range  of  the  A,  B,  C,  and  D  matrices.  There  is  no  direct 
method  of  determing  the  "correct"  values  of  Q  and  R  for  a  particular 
set  of  imposed  control  rate  limits  and  constraints.  The  process  of 
selecting  the  set  of  values  for  Q  and  R  that  results  in  the  desired 
transient  response  is  an  interactive  process.  This  program  and  the 
interactive  graphics  allows  the  control  designer  to  define  the  perfor¬ 
mance  index  and  calculate  the  feedback  gains  for  the  transient  response 
he  desires  in  one  sitting.  It  also  gives  him  the  capability  to 
experiment  with  various  performance  index  weightings  and  view  the 
transient  response  of  the  output  variables  immediately.  He  can  then 
ascertain  the  change  to  the  feedback  gain  matrix 
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As  an  example,  for  the  linear  model  at  PLA  =•  60  (Table  6),  the 
Q  and  R  matrix  values  in  Table  8,  case  1,  were  finally  selected  to 
provide  the  "best"  transient  response.  The  transient  response  as  seen 
on-line  on  the  graphics  unit  is  presented  in  Figure  8,  Case  1.  The  trans¬ 
ient  is  for  a  step  change  in  PLA  from  60  to  63.  The  output  parameters 
displayed  have  been  nond i mens i ona 1 i zed  and  scaled.  In  this  particular 
PLA  operating  range  (CMVT)  it  is  desired  to  hold  the  far.  and  compressor 
speed  and  pressure  ratios  constant,  while  fuel  flow  is  added  to  increase 
turbine  temperature  and  thrust.  To  accomplish  this,  large  Q  matrix 
weightings  were  given  to  the  Nl,  N2, and  WF  errors  with  an  order  of  mag¬ 
nitude  less  weighting  to  the  FN,  TT*iM,  WFAN,  C!VV,and  RCVV  errors. 

Heavy  weights  were  given  to  the  CIVV,  RCVV, and  AJD  rates  of  movements 
in  the  R  matrix.  This  set  of  Q,  and  R  matrix  values  is  one  possible 
performance  index  (PI)  which  gives  one  particular  set  of  feedback  gains. 

As  a  comparison.  Case  2  shows  the  transient  response  for  the  PI  where 
all  the  output  and  control  parameters  were  weighted  equally.  Transient 
response  is  very  sluggish  but  Nl  and  N2  remain  constant.  Case  3  shows 
the  transient  response  for  the  PI  where  very  fast  thrust  response  is 
desired  but  no  other  restraints  placed  on  the  other  parameters.  The 
thrust  response  for  this  case  is  extremely  fast  but  a  very  large  TT4M 
overshoot  occurs  which  is  very  damaging  to  the  turbine.  Nl  and  N2 
move  a  considerable  amount  from  their  respective  CMVT  values.  Although 
optimal  control  techniques  do  not  provide  a  direct  method  to  calculate 
the  "best"  transient  response,  the  interactive  Ricatti  equation  solver 
provides  the  control  designer  with  an  excellent  design  tool  that 
quickly  calculates  a  set  of  feedback  gains  for  a  multivariable  controller. 


Response  for  Various  Performance  Indexes  (PI) 


5.  OPTIMAL  CONTROL  LAW  IMPLEMENTATION  ON  NONLINEAR  DYNAMIC  SIMULATION' 

Before  discussing  the  nonlinear  application  aspects  of  the  optima. 

control  law,  a  summary  is  provided  of  the  steps  taken  to  derive  the 
optimal  control  law. 

a.  Develop  nonlinear  dynamic  model 

b.  Derive  steady  state  control  schedules 

c.  Generate  linear  models  (A,  B,  C,  D  matrices)  at  selected 
power  settings 

d.  Reduce  order  of  linear  models 

e.  Transform  A,  B,  C,  and  D  matrices  to  the  form  of  subsection  3 

f.  Calculate  feedback  gain  matrices  for  small  perturbation  at 
each  selected  power  settings  using  the  reduced  order  models 

A  number  of  problems  arise  in  adapting  the  optimal  control  law 
of  the  last  section  to  the  nonlinear  dynamic  simulation.  These  are: 

a.  Feedback  gains  generated  for  small  transients  won't  work  for 
large  transients 

b.  Control  variable  saturation 

c.  Power  range  of  applicability  for  particular  optimal  control  law 

d.  Transition  from  one  set  of  feedback  gains  to  another  set 

e.  Steady  state  control 

f.  Control  complexity 

The  linear  models  used  for  this  study  were  calculated  about  various 
power  settings  using  small  perturbations.  The  validity  of  a  linear 
model  over  a  range  of  power  settings  has  a  direct  bearing  on  the  range 
of  applicability  of  the  optimal  control  feedback  gains.  At  the  sea 
level  static  condition,  two  linearizations  were  required  to  adequately 
represent  the  VCE  nonlinear  model  over  the  entire  power  range.  Linear 
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models  were  initially  calculated  about  five  different  power  settings 
(PLA  =  25,  35,  50,  60,  83).  Recalling  that  the  VCE  operates  in  two 
modes,  CMVT  above  PLA  =  55  and  non-CMVT  below  PLA  =  55,  an  attempt  was 
made  to  represent  the  entire  power  range  with  just  two  linear  models. 

The  obvious  advantage  of  doing  this  is  the  reduction  in  the  number  of 
feedback  gains  that  must  be  stored.  The  two  linearizations  selected 
were  at  PLA  =  35  and  60.  These  linear  models  are  presented  in  Table  5. 
The  two  linearizations  were  selected  because  their  feedback  gains 
provided  what  was  considered  to  be  adequate  transient  response  over 
their  respective  non-CMVT  and  CMVT  ranges.  Obviously  the  optimal 
feedback  gains  will  not  be  optimal  over  the  entire  operating  range, thus 
resulting  in  a  somewhat  degraded  transient  response.  This  difference 
was  not  quantified. 

Figure  9  shows  a  schematic  of  the  regulator  controller  as  imple¬ 
mented  on  the  nonlinear  transient  simulation.  All  of  the  states,  X, 
are  fed  back  and  an  error  calculated  between  the  operating  value  and 
a  reference  value.  The  reference  values  are  the  steady  state  values 
of  the  control  inputs  and  states  as  shown  in  Figure  5.  The  transformed 
state  matrix,  X,  was  schematically  reduced  to  x  and  i1  to  demonstrate 
that  only  the  four  states  (Nl,  N2,  PT7,  PT6C)  require  measurement. 

The  seven  control  inputs  are  calculated  within  the  controller  so  that, 
assuming  suitable  positioning  mechanisms,  they  do  not  require  measure¬ 
ment. 

The  rate  limited  power  level  angle  (PLAP)  shown  in  Figure  9  was 
Implemented  to  account  for  large  rapid  step  changes  in  PLA.  The  optimal 
feedback  gains  were  calculated  for  small  transients  (PLA  =  60-63). 

For  a  set-point  controller,  the  control  input  is  proportional  to  the 
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Figure  9.  Schematic  of  Regulator  Control 


gain  multiplied  by  the  error  signal.  For  small  PLA  changes,  the 
difference  between  the  state  reference  and  the  operating  values  is 
small  and  the  small  error  signal  multiplied  by  a  given  gain  provides 
suitable  performance.  For  large  PLA  changes,  the  given  gain  multi¬ 
plied  by  a  large  error  signal  causes  large  excursions  in  the  control 
inputs  which  can  cause  temperature  overshoots  and  compressor  stalls. 

PLA  rate  limiting  provides  a  simplified  means  by  which  the  error 
signal  size  is  regulated  by  controlling  rate  of  change  of  the  state 
reference  values.  PLAP  is  allowed  to  vary  at  hO  degrees  per  second. 

This  value  was  found  to  be  the  best  compromise  between  the  small  and 
large  transient  responses.  Other  approaches  to  the  error  magnitude 
problem  include  the  use  of  acceleration  schedules  on  critical  states 
(Reference  7)  and  scaling  the  gains  as  a  function  of  actual  fan  rotor 
speed  divided  by  the  reference  fan  speed  (Reference  6).  All  of  these 
approaches  result  in  the  actual  controller  being  to  some  degree,  degraded 
from  the  optimal . 

Another  aspect  of  the  large  transient  problem  is  the  necessity 
to  transition  from  one  set  of  feedback  gains  to  the  other.  The  CMVT 
breakpoint  occurs  at  PLA  =  55-  One  set  of  feedback  gains  was  derived 
for  the  below  CMVT  range  and  the  other  for  the  CMVT  mode.  A  simple 
approach  to  transitioning  would  be  to  select  a  fan  or  compressor  rotor 
speed  at  which  the  transition  would  occur.  This  approach  was  used  in 
References  5  and  6.  A  result  of  using  this  approach  is  that  discontin¬ 
uities  occur  in  the  output  and  state  variables.  The  approach  used  in 
this  study  was  to  set  up  a  compressor  speed  range  during  which  the  two 
sets  of  feedback  gains  are  gradually  transitioned.  This  logic  is  shown 
schematically  in  Figure  10,  The  method  is  basically  an  averaging  of 
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the  feedback  gains  over  a  transition  range  of  speeds  and  did  an  adequate 
job  of  reducing  the  discontinuity  caused  by  the  transition  from  one  set 
of  feedback  gains  to  another  during  large  transients. 

Control  variable  saturation  presented  an  additional  suboptimal 
factor.  The  development  of  the  optimal  control  law  assumed  that  the 
control  inputs  were  unconstrained.  The  actual  inputs  do  have  physical 
limits,  such  as  minimum  and  maximum  nozzle  areas  or  vane  movements. 

During  a  transient  this  is  not  a  serious  problem  because  the  controller 
will  eventually  bring  the  control  input  off  its  limit.  The  controller 
is  doing  the  best  it  can  by  staying  on  the  limit  as  long  as  the  opti¬ 
mal  request  is  forcing  the  control  in  that  direction. 

When  a  steady  state  condition,  as  determined  by  lack  of  control 
variable  movement,  is  reached  after  a  transient,  control  variable 
saturation  presents  a  fundamental  problem.  The  saturated  control 
variable  is  unable  to  move  to  zero  out  its  own  error,  thus  causing 
the  other  control  variables  to  move  off  their  desired  steady  state 
position,  resulting  in  steady  state  errors.  This  problem  was  dis¬ 
cussed  in  Reference  7-  No  general  solution  was  found  but  a  specific 
solution  was  stated.  Some  of  the  causes  of  this  control  saturation 
include  inaccuracies  in  the  scheduling  and  sensing  of  parameter, 
changes  in  component  efficiences  or  installation  effects  such  as  customer 
bleeds.  These  effects  were  studied  and  are  presented  in  the  next  section. 

Steady  state  control  is  accomplished  using  the  same  feedback  gains 
as  used  for  the  transients.  No  attempt  was  made  to  generate  special 
steady  state  control  logic  such  as  maximum  turbine  inlet  temperature 
limit  or  rotor  speed  schedules.  The  engine  operating  state  is  determined 
by  the  reference  schedules,  installation  effects,  ambient  conditions. 


and  component  eff Iclencles,  Present  day  turbine  engines  require  per¬ 
iodic  adjustments  (trim)  to  their  steady  state  control  schedules  to 
account  for  component  deterioration.  These  trims  usually  affect  only 
the  maximum  (PLA  =  83)  operating  condition  by  moving  the  control 
limits  or  schedules  up  or  down.  For  the  linear  regulator  controller, 
the  reference  schedules  could  be  trimmed  to  accomplish  the  same  effect. 

The  other  approach  would  be  to  allow  the  regulator  controller  to  oper¬ 
ate  with  whatever  steady  state  errors  that  fall  out  and  hope  no  physical 
limits  are  exceeded.  Previous  reports  have  not  adequately  addressed 
this  problem.  The  next  chapter  will  discuss  this  area  for  the  specific 
operating  range  used  in  this  report. 

A  final  problem  area  with  adapting  the  optimal  control  law  concerns 
control  complexity.  By  definition,  linear  regulator  theory  requires 
that  all  the  state  variables  and  control  inputs  must  be  fed  back  to 
have  an  optimal  control.  This  requires  a  good  deal  of  storage  for  the 
feedback  gain  matrices  and  computational  effort.  Various  methods  are 
available  to  reduce  the  number  of  feedback  gain  terms  without  significantly 
reducing  the  system  response  (References  7  and  13). 

The  straight  forward  approach  is  to  compare  the  magnitudes  of  the 
normalized  gains,  eliminating  all  the  gains  less  than  a  preselected  mag¬ 
nitude,  and  reevaluating  response  to  insure  system  response  has  not  been 
significantly  affected.  As  an  example,  for  the  feedback  gain  matrix  of 
Table  9,  CMVT  mode,  core  exhaust  nozzle  area  could  be  calculated  con¬ 
sidering  only  gain  terms  greater  in  absolute  magnitude  than  20.  Instead 
of  11  feedback  terms  only  four  terms  would  be  necessary.  Beattie  in 
Reference  7  reduced  the  number  of  gain  terms  for  his  case  from  50  to  21 
using  this  approach.  A  more  rigorous  technique  was  presented  in 
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Reference  13  which  analytically  evaluates  the  Importance  of  each  gain 
term  on  system  response.  This  was  done  by  evaluating  the  following 
integral  for  each  gain  term  eliminated; 

OO 

PI  =  /  (lY-Y*JTQ[Y-Y*]  +  IU-U*]TR[U~U*])dt 
0 

Where:  Y*  *  full  output  vector 

Y  =  limited  output  vector 
U*  =  full  control  vector 
U  =  limited  control  vector 
Q,  =  output  weighting  matrix 
R  =  control  weighting  matrix 

Although  the  reduction  of  feedback  gain  terms  results  In  a  further 
deviation  from  the  optimal  control  law,  it  is  obvious  that  the  linear 
regulator  theory  and  the  feedback  gain  reduction  methods  provide  a 
systematic  approach  to  determine  the  minimum  number  of  gain  terms 
required  and  their  values.  It  also  provides  a  performance  measure  of 
how  nearly  optimal  the  final  controller  Is.  Conventional  control  tech¬ 
niques  do  not  provide  a  systematic  approach  to  determine  the  cross 
coupled  terms. 


SECTION  IV 


NONLINEAR  MODEL  APPLICATION  RESULTS 

The  previous  chapter  discussed  the  -problem  areas  associated  with 
adapting  the  optima!  control  law  of  Equation  30  to  the  nonlinear  model. 
This  chapter  will  discuss  the  results  from  the  nonlinear  adaptation. 

All  results  are  for  sea  level  static  condition  only,  using  the  two  sets 
of  feedback  gains  given  in  Table  9.  These  two  sets  of  feedback  gains 
were  derived  for  two  unique  performance  indexes  (PI)  that  were  considered 
to  provide  the  "best"  transient  response  while  operating  within  all  the 
constraints.  These  constraints  were  the  same  (turbine  temperature  limits, 
stall  margins  etc.)  that  have  to  be  observed  independent  of  which  design 
method  is  used  to  calculate  the.  feedback  gain  terms.  The  purpose  of  this 
section  is  to  present  results  when  transient  inputs  are  applied  to  the 
nonlinear  model  and  also  to  discuss  how  well  the  control  handles 
"real  world"  transient  and  steady  state  changes  in  the  model,  sensors, 
and  actuators. 

Figure  11  shows  the  transient  response  of  the  VCE  nonlinear  model 
using  the  controller  of  Figure  9.  All  of  the  output,  state, and  control 
variables  are  plotted.  The  transient  is  from  idle  thrust  (PLA  =  20)  to 
maximum  nonaugmented  intermediate  thrust  (PLA  «  83).  Thrust  (FN) 
reaches  S5t  of  its  intermediate  value  in  about  2  seconds  with  no  over¬ 
shoot  or  steady  state  oscillation.  Turbine-in  temperature  (TTAM)  aid 
not  exceed  its  intermediate  value  which  is  critical  to  maintaining 
turbine  life.  The  leveling  off  at  l  second  is  caused  by  the  transition 
from  the  low  feedback  gain  matrix  to  the  high  gain  matrix.  Fan  airflow 
(WFAN) ,  which  is  important  for  spillage  drag  as  explained  earlier, 
reaches  steady  state  after  only  1.25  seconds.  Because  thrust  is  a  com¬ 
bination  of  airflow  and  pressure,  this  says  that  the  first  half  of  the 


transient  thrust  Increase  is  contributed  almost  entirely  by  airflow 
and  fan  duct  pressure  (PTbc)  Increase.  Core  exhaust  pressure  (PT7) 
remains  almost  constant  the  first  second  and  then  increases  rapidly 
for  the  next  second.  As  a  comparison,  Figure  16  (base  case  displayed  later) 
shows  the  same  transient  for  the  original  FI 00  nonlinear  model.  Thrust  levels 
off  at  93%  at  about  3-5  seconds  but  does  not  reach  full  thrust  until 
about  5.75  seconds.  Maximum  airflow  is  not  reached  until  about  2.75 
seconds.  The  stall  margins  (SMAF  3nd  SMHC)  are  approximately  the 
same  for  both  models.  V/einberg  (Reference  6)  using  the  original  FI 00 
nonlinear  model  and  a  multivariable  control  attained  maximum  thrust  and 
airflow  In  about  3  seconds.  Beattie  (Reference  1)  using  a  VCE  model  and 
conventional  cjntrol  design  showed  thrust  and  airflow  transients  of  5 
seconds  to  attain  maximum  values.  Using  the  same  model  but  with  a 
multivariable  control,  Beattie  (Reference  7)  showed  a  thrust  response  of 
h  seconds  for  the  idle  to  intermediate  transient,  but  for  a  small  tran¬ 
sient  similar  to  Figure  1 3 >  the  transient  response  occurs  in  0.25  seconds . 

Figure  12  presents  a  deceleration  from  intermediate  thrust  to 
idle  thrust.  The  time  required  to  reduce  thrust  is  an  important  con¬ 
trol  design  consideration  for  air  combat  manuevers.  idle  thrust  is 
reached  after  about  2.5  seconds  with  no  problem  with  stall  margins. 

The  only  significant  problem  was  caused  by  the  fuel  flow  (WF)  going 
too  low  and  possibly  causing  a  flameout.  This  could  be  overcome  by 
putting  a  minimum  fuel  flow  limit  on  the  control.  This  would  not  affect 
the  thrust  response.  Maximum  and  minimum  limits  were  placed  on  some 
of  the  other  control  acutators.  The  dashed  lines  in  Figure  5  show  these 
i imits. 

Figure  13  shows  a  small  transient  in  the  CMVT  range.  In  this  range 
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Figure  13-  (concluded) 


rotor  speeds  remain  essentially  constant,  fuel  flow  Is  Increased  to 
raise  turbine  temperature,  the  turbine  nozzles  (HTVPOS  and  FTVPOS)  are 
positioned  to  maintain  compressor  and  fan  operating  points, and  core 
exhause  nozzle  (AJC)  is  closed  down  to  build  up  PT7  and  increase  thrust. 
Intermediate  thrust  is  reached  in  about  one  second. 

Figures  11,  12, and  13  show  clearly  that  the  multivariable  con¬ 
troller  described  in  this  report  is  capable  of  controlling  the  VCE 
nonlinear  model  during  large  or  small  transients  and  steady  state,  as 
well  or  better  than  the  less  complicated  original  F100  engine.  The 
multivariable  controller  was  designed  for  a  single  operating  condition 
and  adapted  to  a  nonlinear  model  from  which  the  linear  models  were 
derived.  In  an  actual  application  the  controller  would  have  to  operate 
at  many  different  flight  conditions  on  an  engine. that  might  behave 
differently  than  the  nonlinear  model  and  have  mlsrigged  actuators  or 
miscal ibrated  sensors.  To  study  how  these  factors  are  handled  by  the 
multivariable  controller,  transient  responses  (PLA  «*  20  to  PLA  =  83) 
were  run  for  various  cases  of  implanted  faults.  The  steady  state 
results  are  tabulated  in  Table  10  and  discussed  below.  The  affect  on 
transient  response  for  all  the  cases  was  insignificant.  For  comparison 
two  of  the  cases  were  run  on  the  original  F100  nonlinear  model. 

All  engine  controls  require  some  sensor  inputs  so  the  controller 
can  properly  position  the  control  actuators.  Of  importance  to  the  con¬ 
trol  designer  is  the  impact  of  sensor  errors  on  critical  output  para¬ 
meter  response.  The  multivariable  controller  presented  in  this  repori 
uses  four  sensed  values  (Nl,  N2,  PT7,  PT6c) .  These  four  sensed  values 
along  with  the  seven  calculated  control  actuator  values  are  fed  back  to 
calculate  the  rate  of  change  of  the  control  actuators.  To  simulate  the 


effect  of  a  sensor  error  the  sense  value  of  N2  and  PT7» going  Into  the 
control ler,  were  multiplied  by  0.97  one  at  a  time.  The  multivariable 
controller  now  has  a  3%  error  in  one  out  of  eleven  values  that  deter¬ 
mine  rates  of  change  of  the  seven  control  actuators.  If  the  actuators 
are  adjusted  to  zero  o=t  the  error  in  this  one  sensor  then  errors 
would  appear  in  the  other  ten  feedback  terms.  If  the  gain  term  for  the 
sensor  with  the  error  is  small  compared  to  the  sum  of  the  other  10  gain 
terms  the  effect  on  controller  output  will  be  minimal.  For  the  two 
sensors  selected  this  was  the  case.  Table  10  shows  minimal  changes  in 
the  output  parameters  for  N2  and  PT7  sensors  errors  of  3%- 

Another  consideration  for  the  control  designer  is  actuator  posi¬ 
tioning  errors.  These  errors  tend  to  have  a  greater  impact  on  the 
output  parameters  becruss  the  feedback  gain  matrix  weights  the  actuator 
terms  more  than  the  state  sensor  terms.  In  fact,  one  of  the  basic  assump¬ 
tions  made  in  the  development  of  the  multivariable  control  was  that  the 
seven  control  actuator  positions  did  not  have  to  be  measured.  The 
calculated  position  is  fed  back  and  the  actual  positioning  is  assumed 
to  be  correct.  To  determine  the  affect  on  the  output  parameters  for 
actuator  positioning  errors,  the  nonlinear  model  was  run  with  various 
positioning  errors  introduced.  The  actuator  errors  simulated  were  fuel 
flow  (WF),  fan  duct  nozzle  area  (AJD) ,  high  turbine  nozzle  .’rea  (HTVPOS)  t 
and  compressor  vane  position  (RCVV).  The  errors  were  simulated  by  sub¬ 
tracting  the  amount  shown  in  Table  10  from  the  calculated  position  and 
feeding  this  back  to  the  controller.  WF  and  AJD  errors  had  the  largest 
effect  on  the  output  parameters.  A  300  pound-per-hour  error  in  WF  is 
equal  to  about  33  at  Intermediate  power  setting  and  the  resulting  change 
in  thrust  is  not  unreasonably  large.  The  original  F100  control  schedules  a 


turbine  temperature  and  adjusts  WF  to  hold  tnis  temperature.  Therefore 
a  fuel  flow  error  would  have  no  affect.  But  the  temperature  measure¬ 
ment  accuracy  is  +0*5%  which  can  cause  a  variation  in  FN  of  about  1.5%. 
The  RCVV  error  was  handled  with  little  difficulty  by  the  multivariable 
controller.  The  full  transient  response  for  this  case  is  presented  in 
Figure  U  with  the  base  transient  for  PLA  -  20  to  83.  For  comparison 
Figure  If-  shows  the  effect  of  this  same  error  on  the  transient  and 
steady  state  response  of  the  original  F100  model.  As  shown  in  Table  10 
this  error  has  a  very  large  effect  on  the  output  parameters. 

A  basic  assumption  of  the  multivariable  control  is  that  the  linear 
model  used  to  derive  the  feedback  gains  accurately  represents  the 
actual  system.  Some  of  the  mechanisms  that  can  change  the  model  are 
component  efficiency  deteriorations,  changes  in  airflow  pumping  char¬ 
acteristics,  pressure  loss  changes,  or  externa!  extractions  of  compressor 
bleed  or  horsepower.  Current  turbine  engine  controls  require  some  sort 
of  retrimming  of  one  or  more  control  schedules  when  component  changes 
occur.  External  extractions  usually  are  not  accounted  for  and  result 
in  changes  to  the  output  parameters.  Two  of  these  mechanisms  were  sim¬ 
ulated  on  the  VCE  model  and  or.e  of  them  on  the  original  FI 00  model  with 
the  results  presented  in  Table  10.  In  one  case  high  turbine  efficiency 
was  reduced  by  2%.  The  VCE  multivariable  control  handled  this  change 
to  the  model  wit!  only  minor  changes  to  the  output  parameters  whereas 
the  original  F100  control  allowed  FN  to  change  by  **.3%.  The  transient 
responses  for  a  high  turbine  efficiency  change  for  both  models  is  shown 
in  Figures  15  and  17.  The  other  mecha>.:«m  studied  was  compressor  bleed 
air  extraction  of  1  lb. /sec.  Again  the  multivariable  control  was  able 
to  handle  this  with  considerably  less  effect  on  FN  than  the  original 
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FI 00  control  (0.9%  change  vs.  2.2%). 

As  a  final  item  Figure  18  shows  the  response  of  the  output  para¬ 
meters  for  a  reduced  feedback  gain  matrix.  This  reduction  was  discussed 
in  Subsection  1 1 1-5  and  was  carried  out  to  show  to  what  extent  the  optimal 
control  law  gains  could  be  reduced  and  still  maintain  the  same  transient 
response.  The  original  feedback  gain  matrix  had  77  terms,  while  the 
reduced  matrix  had  only  33  with  no  change  in  the  output  parameters. 

The  method  used  to  reduce  the  number  of  terms  was  to  calculate  the 
percentage  contribution  of  each  term  on  the  seven  different  actuators 
and  retain  the  largest  terms  that  added  up  to  a  preselected  percentage. 
This  method  required  little  trial  and  error  but  further  reductions  could 
probably  be  obtained  using  the  method  detailed  in  Reference  13. 
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Reduced  Feedback  Gain  Matrix  Effect 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  demonstrated  the  usefulness  of  linear  quadratic 
regulator  theory  as  a  design  tool  in  developing  a  controller  for  a 
multi-input  multi-output  system.  Without  such  a  tool  it  would  be  very 
difficult  and  time  consuming  to  calculate  all  of  the  feedback  gain 
terms,  accounting  for  all  the  cross  coupling  effects  neces.  ry  to  take 
full  advantage  of  the  additional  control  input  variables  provided  by 
a  variable  cycle  engine  (VCE)  configuration.  This  is  nrt  to  say  that 
a  controller  developed  using  classical  control  techniques  could  not 
provide  the  same  transient  and  steady  state  performance,  but  that  the 
time,  effort  and  intuition  required  can  be  significantly  reduced  using 
the  approach  described  in  this  report. 

The  interactive  graphics  program  developed  in  this  report  solves 
the  RicattI  equation  for  a  user  chosen  performance  index,  calculates 
the  feedback  gain  matrix  and  displays  the  transient  response  of  the 
output,  state  and  input  variables  on-line.  This  greatly  reduces  the 
time  and  guessing  required  to  develop  values  for  multivariable  gain 
terms  by  providing  a  systematic  method  for  calculating  all  of  the 
cross-coupling  effects.  In  addition,  the  control  designer  has  a  good 
idea  of  what  the  actual  transient  response  of  the  engine  will  be. 

The  controller  developed  in  this  report  is  capable  of  controlling 
the  VCE  nonlinear  model  for  both  transient  and  steady  state  at  the  sea 

level  operating  condition,  for  both  large  and  small  step  changes  in  the 
requested  thrust  level  and  within  the  same  constraints  imposed  on  the 
original  F100  model.  The  transient  response  times  were  considerably 
better  than  the  original  F100  model  owing  to  both  the  multivariable 
control  concept  and  the  inherent  capabilities  of  the  VCE  concept.  The 
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transient  responses  were  similar  to  published  results  for  similar 
engine  cycles.  The  steady  state  variations  In  the  output  due  to 
"real"  world  effects  were  also  handled  considerably  better  by  the 
multivariable  controller  than  by  the  original  F100  model.  This  was 
due  to  the  different  set  of  control  variables  used  by  the  multivariable 
controller  and  the  inherent  capability  of  the  multivariable  controller 
to  compensate  for  variations  in  the  system. 

The  controller  developed  in  this  study  was  suitable  only  for  the 
sea  level  static  condition.  Further  work  should  be  directed  toward 
the  calculation  of  gain  matrices  that  would  cover  the  entire  flight 
envelope.  This  would  require  improved  gain  transitioning  algorithms 
and  steady  state  control  schedules.  Also,  further  studies  into  the 
many  "real  world"  mechanisms  that  customarily  affect  turbine  engine 
o  eration  should  be  performed.  A  method  will  also  have  to  be  developed 
which  allows  optimization  of  the  feedback  gains  on  the  nonlinear  model 
and  the  actual  engine.  That  is,  if  the  feedback  gain  matrix  calculated 
from  linear  models  does  not  control  the  actual  engine  properly,  a 
method  must  be  available  to  determine  which  terms  of  the  gain  matrix 
have  to  be  changed  and  by  how  much.  A  future  study  effort  should  also 
determine  if  a  turbine  temperature  sensor  is  necessary  to  safely  control 
a  turbine  engine.  Sensor  and  actuator  dynamics  should  also  be  simulated 
to  better  represent  the  actual  engine. 
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APPENDIX  A 


DERIVATION  OF  NOZZLE  POSITIONING 
FOR  VARIATIONS  IN  TURBINE  INLET  TEMPERATURE 

Two  basic  equations  will  be  used  to  derive  relationship.'  for  the 
three  areas  (A4 ,  A45,  and  A8,  equivalent  to  HTVPOS,  FTVPOS,  and  A-JC)  as 
T4  (1T4M)  is  varied  to  modulate  thrust  while  a  constant  match  of  the 
fan  and  compressor  is  maintained.  Airflow  (WG4) ,  pressure  (PT4)  and  the 
work  (Ah)  across  the  turbines  will  be  nearly  constant  because  of  the 
constant  match  of  the  fan  and  compressor. 


T4 

— Sf - 

- v - 

- ST 

PT4  «♦ 

A4 

A45 

A8 

WG4-> 

-A _ 

_ A _ 

_ A _ 

Starting  at  Station  k  (high  pressure  turbine)  and  assuming  for 
simplicity  the  turbine  efficiencies  (n)  remain  constant,  the  first  two 
equations  are: 

-  WG4  (T4)15  /  (PT4  Ak)  =  K1  (A"D 

Ah  =  n  Ah..  «  n  c  n  (1  -  (PTi*5  /  PT4)*Y~^/y)  =  K2  (A-2) 

I O  P 

»  n  C  TMl  ■  T45  /  T4)  =  K2 

The  relationship  for  A4  is  obtained  directly  from  Equation  A-l; 

A4  »  K3  (T4)*5  where  K3  =  WG4  /(PT4  K1)  (A"3) 

The  relationships  between  T4  and  Pt45  and  T45  are  obtained  by 
rearranging  Equation  A-2; 

PJks  /  pt4  -  (l  -  K2  /  (n  r  t4))y/(y_1) 
or  p 

PT^*5  -  PT4  (1  -  Kk  /ti»)y/^y"^  where  Kk  *  K2  /(nCp)  (A-  • 
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(A-5) 


Rearranging  and  substituting  for  PT45  /  PT4; 

ns  -  t4  -  K2  /  Cn  c  )  »  t4  -  k4 

p 

Continu:ng  to  station  45  (fan  drive  turbine)  and  using  the  same 
basic  equation; 

r45  =  WG45  (T45)^  /  (PT45  A45)  =  K5 
Ai.  =  n  Ahfd  =  n  Cp  T45  (1-  (PT5  /  PT45) (y~1)/y)  =  K6 
»  ti  Cp  T45  (1  -  T5  /  T45)  =  K6 

Rearranging  EquationA*6and  substi tuting  ;or  PT45  and  T45  from  Equations 
A-4  and  A-5,  A45  is  obtained; 

A45  -  WG45  (T45)**  /  (PT45  K5)  =(K7  (T4-K4)^  T4y/(y“,) )/(PT4  (t4-k4)y/(y_1)) 

where;  K7  =  WG45  /  K5  (A-8) 

PT5  is  obtained  from  Equat ionA-7 by  subst i tut ing  in  Equations  A-4  and  A~5 
PT5  /  PT45  =((1  -  K6  /  (n  Cp  T45))y/(y"1; 

or 

PT5  =  (1  -  K4  /  T4)y/(y_1)((1-  k8)  /  (T4  -  K4))y/(y“^  (A-9) 

where  K8  =  K6  /  (n  Cp) 

T5  =  T45  -  K6  /  (n  Cp)  =  T4-  K9  where  K9  =  K4  +  K8  (A- 10) 

The  relationship  for  the  core  nozzle  area  A8  is  derived  assuming  a 

choked  nozzle  and  a  constant  pressure  drop  across  the  tailpipe. 

r8  -  (WG8  T8li)  /  (PT8  A8)  =  K10  WG8  is  constant  (A_n) 

PT8  =  PT5  (1  -  Kll) 

T8  =  T5 

rearranging; 

A8  »  K12  T5 *  /  PT5  where  K12  =  WG8  /  (K10  (1-K11)) 

A8  «  (K12  (T4  -  K9)^)  /  (PT4  (1-k4/T4)y/^y"1  ^  (1-k8/<T4-K4j)y/^y"^) 

(A- 12) 

Assuming  some  typical  values,  for  the  consents,  the  three  areas 
can  be  calculated  for  various  values  of  T4. 


(A-6) 

(A-7) 
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WG4  -  117.2 

WG45 

«  140.9 

PT4  «  354,7 

K5 

-  67.1 

K1  «  18.2 

«  .903 

tv4  =  .888 

V5 

-  .29 

C  4  «»  .30 

F 

K6 

-  97.6 

K2  »  227. 

WG8 

*  141.7 

Y  »  1 .345 

K10 

-  132.7 

Kll 

=  .038 

Using  these  values,  Equations  A-3,A-8  and  A- 1 2  were  evaluated  for 
various  values  of  T4  and  plotted  in  Figure  A1  with  the  actual  areas  from 
the  nonlinear  simulation.  The  areas  are  relative  to  their  intermedia's 
level  values. 
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Figure  A-l.  Variation  in  Areas  with  Turbine  Temperature 
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